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Abstract 

We consider a Callan-Symanzik and a Wilson renormalization group (RG) 
approach to the infrared problem for interacting fermions in one dimension 
with backscattering. We compute the third order (two-loop) approximation 
of the beta function using both methods and compare it with the well known 
multiplicative Gell-Mann Low approach. We point out a previously unnoticed 
strong instability of the third order fixed point with respect to an arbitrary 
dimensionless parameter, which suggests a RG flow toward a strong coupling 
phase. 
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I. INTRODUCTION 



The problem of the one-dimensional Fermi gas model of a metallic conductor, in the 
low energy approximation, has been approached using three methods: conventional many- 
body techniques^ and, mainly, the bosonizationlHl and the renormalization group (RG) 
methodlHBra. In this paper we will be concerned with the latter approach. A formu- 
lation of the Gell-Mann Low multiplicative RG for this problem was introduced in Ref. 
|5|. The model considered was the g-ological model, which describes a weakly interacting 
one-dimensional fermion system with Tomonaga-type (g 2 , 9i) and backscattering (gi) in- 
teractions. Phonons are neglected. That method provided a satisfactory understanding of 
the infrared behavior in the case of a weak repulsive (effective) interaction. A short list of 
the most relevant results in this case may be the following (for extensive reviews see e.g. 
Refs. 0,^): the RG flows toward the Luttinger liquid^ fixed point00; ii) there is a line 
of nontrivial fixed points; iii) in the infrared limit the system is not asymptotically free, 
as in the Fermi liquid case, but is described by anomalous exponents. These results were 
recovered and rigorously proved also in the case of periodic potential using the Wilson rgO. 

Things change considerably if we consider a weak attractive interaction. Since in this 
case there is not a second order (one loop) finite fixed point, in Ref. § the computation of the 
beta function was carried to third order (two loops). It was found a 0(1) third-order fixed 
point. This result, if reliable, would be of extreme physical interest because it would signal 
a behavior completely different from the Luttinger liquid paradigm. One should expect the 
opening of a gap in the dispersion relations, while Luttinger spectrum is gapless, and an 
exponential decay of the correlation functions, while in the Luttinger case there is only a 
power low decay with increasing distance. The problem is, of course, how seriously one 
should take the very existence of an attractive perturbative fixed point on the basis of the 
third order result. The computation of the the fourth order (three-loop) approximation of 
the beta function was discussed in Refs. pOH22|. A smaller but still 0(1) fixed point was found. 
Moreover only the first two terms of the beta function are universal. The computation of 
the third term is useful provided there is some evidence of a perturbatively tractable phase 
interacting attractively. In this case a precise determination of the renormalized couplings 
would be important to compute the response functions. 

It it useful to make a comparison with the results obtained with the bosonization method. 
With bosonization it is meant the bosonic representation of fermion field operatorsll'il'ii. 
This method is in some sense the inverse of the one used to solve exactly the Luttinger 
model@, where bosonic degrees of freedom are expressed in terms of fermionic operators. 
Probably the most important result of the bosonization is the exact solution of the model 
with backscattering! {gi and g 2 terms, see below) in the particular case where g\ = — §tt- 
Actually the decoupling between charge and spin degrees of freedom, crucial for the exact 
solution, is open to questions!!. Moreover there are problems in the limiting procedure 
employed^ and the ladder operators restoring the correct occupation numbersB are not 
discussed. (A version of the bosonization free from this problems has been proposedil. It 
should be noted that this version does not deal with the crucial backscattering interaction 
term: in Ref. ^ only the Luttinger model is considered). Anyway, taking for granted the 
Luther-Emery solution!, the RG method should fill the missing information for value of 
gi in the neighborhood of the exact solution. From the bosonized representation of the 
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interaction it is not difficult to derive the third order scaling equationsiil and the response 
functions calculated in Ref. 20 are in good agreement with the results of Ref. ||. 

From these considerations one may be tempted to give a heuristic meaning to the large 
but finite fixed point. In this paper we want to show that this is not the case. The main 
point is that even the sign of the third order fixed point depends on small variations of a 
parameter 7 whose value can be arbitrarily chosen, provided 7 > 1. We will show this both 
using the Gell-Mann Low (GML) and the Wilson RG. 

The paper is organized as follows. In section II we briefly review the multiplicative GML 
approach. We explain why it is useful to check the results of this approach using other 
methods. Recasting the multiplicative procedure into discrete steps, instead of consider- 
ing the usual Lie equation, we reach our main conclusion. In section III we formulate a 
Callan-Symanzik (CS) approach to the problem and compute the beta function in two-loop 
approximation. The same computation is proposed in section IV employing the Wilson RG 
in the multiscale formulationci'B. Finally in section V we come to the conclusions. 



II. THE GELL-MANN LOW APPROACH 

We briefly recall the GML multiplicative RG for one dimensional interacting fermions. 
We will follow closely Refs. [5],|6], with the only difference that we find it convenient to adopt 
a Euclidean formalism. We consider the g-ological model, defined as follows. The kinetic 
term is taken linear around the Fermi surface defined by the two points — hp and kp\ 

where ip k w a are creation and annihilation operators for right moving {u = 1) and left moving 
(u = — 1) fermions with momentum k and spin a [a = ±1/2). We choose units such that 
Vp—1 (vp is the velocity at the Fermi surface). The ultraviolet (u.v.) stability is imposed 
by bandwidth cutoffs: the momenta are restricted to the intervals (ukp — k uv ,uk F + k uv ) 
for i/)^ ff . We define E = 2k uv . The interaction Hamiltonian is 

Hint= 2£ (fl , l|| < W + 9l±^a-a')'iPk 1 ,uj,a' l l ! k2 -uJ,a>i'k2+2k F +p,u} ) a'i'k 1 -2k F -p-uj,a 

k,p,u),a,a' 

+ £ fencer' + 52X^,-«T')^^, <r V'£,-w,^V'j^H^,-w,^^fc 1 - 1 ,^ f «T 

/e,p,oxcr,cr / 
k,p,u),a,a' 

L is the length of the line. The umklapp interaction term (#3) is neglected since it is 
important only in the half-filled band case, which will be excluded. Since gn = —g2\\ it is 
always possible to take g 2 ± = g 2 \\ = #2, reducing the independent couplings to <7i||, gi±, g 2 - 
For the sake of simplicity it is possible to neglect, at least as a first approximation, g 4 : we 
know from the Mattis modefl that g^ does not change the essence of the problem. 
In the Euclidean formalism the free propagator in momentum space is given by 
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G u (k) 



-ik + uki ' 



(2) 



where ko is the energy, k\ the momentum (measured from the Fermi surface), k = (ko, k\) and 
iv — 1 (—1) for right (left) moving fermions. The renormalization procedure is a prescription 
that defines new couplings for a theory with a lowered u.v. cutoff Eq. In the limit Eq — > 
we obtain the renormalized couplings. If Gf is the interacting propagator, the d function is 
defined by the relation 

ki ko 

v kuv Eq / 

The multiplicative constants z and Zi (i — 1 ||, 1 _L, 2) that relate d and the adimensional 
vertex functions Tj for different values of the cutoff are definite by: 



fh k A (E' \,fh k N 



'K 



k' Ei 



Eq 



kj_ 

k U v Eq 



9i 



'K 

Eq 



(3) 



where E' Q < Eq is the lowered cutoff, g and g' denote respectively the old and the new 
couplings. The invariant couplings gf are defined by: 



9? (^,9 



The gf are invariant in the sense that 

E 



( E 



9iZ [—,g)zi(— : g 



( E 



9, 



R 



M' 9 . 



\E 



3 >(f a - g 



\Eq ' 



(4) 



The couplings g\ for the theory with u.v. cutoff E' are defined by 



9i = 9 



Ft I ^0 
\Eq 



,9 



(5) 



A differential equation for gf is readily derived and is the standard Lie equation: 
where x = E' /E . We are interested in the scaling limit x — > O.The two-loop result is 



dg[ 



R 



dx 

R 

1± 



dg R 



1 

x 
1 



-n R2 
9i± 

7T 



1 



2tt 2 



n R n R2 

9i\\9i± 



dx x 
dg R = 1 

x 



1 .R,R 1 



dx 



-9n9i± + ^M\9i± + 9i±) 

R R7 



I I 



(6) 
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For spin independent interaction (g^ = g±± = g{) the nontrivial fixed point is found for 
gt = -2n. 

We now want to recover this result iterating by discrete steps the procedure that defines 
the new couplings when the cutoff is lowered: we aim to study the dependence on the scaling 
parameter. Let 7 > 1. In proper units we put E = 7 and g^ = gi(E ) for i — 1 ||, 1 _L, 2. 
gi-i is defined as (see Eq. (||)) 

where g^o) — g{Eo) = g. The procedure is iterated in the following way: for n < we define 

From Eqs. (^) e (^) we have that the gi jTl for n — —1, —2 . . . are the couplings corresponding 
to the cutoff sequence {7™}: 

In the limit n — > —00 we get the renormalized couplings. We have: 

9l±,n , 

0i||,n-i = 9i\\,n In 7 

71 

+^9i\\,n9i±, n (^ln 2 7 - -ln 7 

gi±,n-l = 9l±,n ~ -gi±,ngi\\,n^1 
7T 



^ (gfi n gi±,n + gl^n) ( ln2 7 - \ In 7 



#2^-1 = #2,n - 2-Su,„ln7 



2tt 2 



9i\\,n9i±, n ln 7 - o ln7 • (7) 



It is easily checked that in the limit 7 — > 1 + Eqs. (|6|) are recovered. In general the fixed 
point depends on 7 (let's remember that it is a third order fixed point). For 7 7^ ^/e, if 
01H = 9iJ- = 9i, we have: 

When 7 — ► 1 we obtain the previous result 5^ = — 2tt. The Lie equations (||) should not be 
fundamental and we find no reason to use the continuous RG instead of its discrete version. 
The dependence on 7 will be discussed in section V. 

As a final comment on this method we note that Eqs. (|3]) rely on neglecting small 
contributes that would not allow one to set multiplicative relations where the z factors do 
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not depend on the external momenta. For example in the case of the one-loop approximation 
of the four-point vertex function, proportional to 

(k is the external energy, see figure [[]), the second term is neglected. That is to say that the 
vertex functions can be divided in scaling and not scaling terms. The first ones are taken 
into account while the second are not discussed in |5]|| For this reason we find useful to 
check Eqs. (|J) using other methods. 



III. THE CALLAN SYMANZIK APPROACH 

Within the framework of the multiplicative RG, it is not difficult to formulate a Callan- 
Symanzik approach for our problem. We follow a common procedure: first we renormalize 
the theory in the u.v. with a fixed (renormalized ) i.r. cutoff m, then we will compute the 
beta function and study the i.r. behavior for m —>■ 0. This approach, devised for a Field 
Theory, in our case may be considered unnecessary. Nevertheless we consider it a way to 
support the GML result. 

The i.r. regularized free propagator is defined by inserting a bare mass tuq in the prop- 
agator (g): 

G ( k m 2 ) - i h±^l 

where k 2 = k 2 , + h\ and again uj = 1 (—1) for right (left) moving fermions. We know that 
the Luttinger model with a local interaction is not renormalizable in the u.v.0 (this is also 
seen from the exact solution^). In order to impose the u.v. stability we choose a nonlocal 
interaction whose strength decreases with increasing distance. The interaction Hamiltonian 
of the model is 

E / d2% d2 y ^t,^t-^' v n ( x ~ v^y^Ax-^a 



J2 J d 2 x ^yip+^ip+^V^ix - y)^ y ^Ax-w 

<jj.cr,cr' 

+ J2 J d2xd2 y^At^^ iV ^ x -y)%^,Ax,u,a, 



where ip^u a are the fermion field operators in coordinate space. The potentials V{ may be 
chosen for instance as follows: 

V t (x) = ^pe- p ^5(x ), i = l ||,1 ±,2, 4, (9) 
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where p > is fixed; xo and x\ are the time and space coordinates. In momentum space the 
model is the same as ([!]) with the only difference that the bandwidth cutoffs are replaced 
by the nonlocal couplings 



P 2 



where k\ is the exchanged momentum in the given interaction vertex and p is introduced in 
Eq. (H). In the limit p — > oo we recover the local couplings of Eq. ([[]). 

Fermion loops are logarithmically divergent. The theory is regularized introducing a 
cutoff A by the means of the standard Schwinger parametrization: 



k 2 + rrin 



In order to renormalize the theory we find it convenient to follow the scheme for the local 
(p = oo) case, even if when p is finite we make more subtractions than strictly necessary. It 
is a simple exercise of standard power counting to find the superficial degree of divergence 
D for the n-point vertex functions in the local case: 

D{T n ) = 2--. 

We renormalize the couplings (gi — > g R ) , the mass (m — > m) and the wave function (ip — > 
il) ±R ). The multiplicative constant Z is formally introduced by the relation = Z 1 ^ 2 ip ±R . 
Let be the renormalized proper n-point vertex functions. The relation between bare and 
renormalized vertex functions is: 

r n (q, m , g, A) = Z- n ' 2 Y R ; (q, m, g R ), (10) 

where q denotes the n — 1 independent external momenta of T n and g = {gi\\, gix, #2> 94}- 
Of course the T n are functions of the spin and to indices attached to the external fields. For 
simplicity we have not indicated this explicitly in Eq. flTU|). For the four-point functions the 
different possible cases are labeled by a single index i — 1 ||, 1 J_, 2, 4. It should be noted 
that while r 4 j do not depend on the value of u, T 2 does, so we need an to label for this 
vertex function. It proves useful to introduce the reduced two-point vertex function T 2 {k): 

f 2 (k) = (ik + uk l )r 2iU> (k). (11) 

In the local case T 2 (/c) does not depend on u. In the non local case the non vanishing terms 
in the infrared limit will be u independent, so we will neglect the u dependence of T2(k). 
The normalization conditions, which define g R , m and the finite part (zero-loop term) of Z, 
are: 

r?i(0)=3? (12a) 
ff(0)= m 2 (12b) 

' 9 -^2(k) u =l. (12c) 



2k dk t 







fe=(0,0) 



7 



The CS equations are derived considering insertions of operators related to the derivatives 
of the vertex functions respect to the i.r. cutoff m . To this end we introduce the operator 
O: 

d 2 x O w 



0(*) = £/ 



2n xq — z — iuj(xx — z\) 



The corresponding source term in the action has the form / d 2 xv (x)O(x) with [v] = 2. In 
momentum space the operator O is 



O(q) = E / 

i . 1 /-r 



d 2 k 



(2tt) 2 [i(k + q ) + Lu(h + qi )] 



where ipk,u),a are the field operators in momentum space, and q is the external momentum 
of the inserted O operators. When O is inserted in a vertex function T n , the value of D for 
r n ,o is 

D(T n , ) 



n 



2-- + ([0]-2) 



n 
2' 



The previous relation means that no new u.v. divergences appear due to O insertions (we 
do not consider r .o We remember that the vertex functions with insertions are defined as 
usual by the Legendre transformation on the field source only and not on the source of the 
inserted operators). Since O does not introduce new divergences, we have O = ZO . The 
insertion of s operator O in a vertex function T n will be denoted with r njS . In analogy with 
Eq. (|TTD we define r2, s (<Z, k) and T^siq, k), where k denotes the s external momenta of the 
s inserted O operators. Equation flTU| ) generalizes into 



r n , 8 (q, k, m , g, A) = Z- n ' 2 Z s T* s {<i, k, m, g R ) 
It is easily deduced that 

d 



r„( q ) = r n>1 ( q ,o), 



(13) 



(14) 



where the insertion of O in the r.h.s. is made at zero momentum, as indicated. From (|13|). 
(|TJ) and (|n|) we have 



l?i(0) 



(15) 



m 



Zm 



o- 



where the second relation follows from Tf (0) = ZT2(0). From Eq. (14) we have: 



d 

m— r n (q, m ,g,A) 



dm 2 . 
= m—— 

g ,A am 



g,A 



r n ,i(q, 0,m ,g, A). 



(16) 



From equations ( fL6|) and ([IB]) , with the definitions 

7i 

(3i= m 



1 dZ 

di 
dgf 



Z dm 



dm 



g,A 



S 



(ji, defined in the previous equations should not be confused with the RG rescaling factor 
7) we obtain: 



m 



_d_ 

dm 







Z m 



dm?. 



dm 



r«(q,0,m,^). 



go,A/ 



(17) 



It is easy to eliminate any reference to the bare theory. From Eq. (|T7| ) written for n = 2 we 
have: 



/ dml 
Z m 



dm 



TUq,0,m,g H ) 



go, A/ 



From Eqs. ( |1 2b|) and the first of (|T5D we conclude 

1 dm?. 



Z m- 



dm 



go, Ay 



(2-7lK, 



so that Eq. (O) can be written as 



( m A + X> (9 «) 



d 



71 



dgl 



„ - 7l jr^(q,m,^) 



(2- 7l )m 2 r5(q,0,m,^). 
The generalization of Eq. (|17|) to the case of s insertions is immediate: 



(9m 



-^ + s)Ti)r5 a (q,k,m^ fl ) 



(2- 7l )m 2 r^ +1 (q,k,0,m^ K ) 



Having set the general definitions and relations of the CS approach we can proceed. We 
will limit ourselves to the computation of the beta function, which is our problem. The 
normalization conditions ( |T2| ) fix the zero-loop terms in the loop-wise expansion of gi, mo 
and Z: 

gf= gf ] 

2 2(0) 

m = m 
Z®= 1. 

One-loop calculations are easily done. Of course m^ 1 ' = and Z^ = 0. This implies that 
up to one loop 71 = (that is 71 = 0(g 2 )). It is convenient to write down the results for 
the couplings in terms of g 4 , gi±, g 2 and g = g 2 — gn. We find: 
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o (1) 



-[lnm-(ln2 + lnp-l)]^ 2 R 

7T 



+ - 



7T 



«?>= 



In m — In 2 + In p — 



1 



— [lnm - (In 2 + lnp- 1)] gfl 
2n 



1 

2^ 



In m - In A + - ( 1 + C + In 2) 
2 



2tt 



a (1) - - — a R2 



(18) 



where C is the Euler constant. In Eq. (jig) we note the presence of p. However simply on the 
basis of dimensional analysis we can exclude that p will appear in the final result. Two-loop 
calculations are tedious and we will omit the details. We calculate only the singular terms 
in m since we do not plan to go beyond the two-loop approximation. We report the results 
for Z and gi±: 



Z {2) = In m (g R ± 2 + gf + 2g^ - 2g%g« 



(19) 



(2) 

9x,± 



+ 



2tt 2 



2vr 2 



lnm — 



In 2 m — 2 lnm (In 2 + lnp — 1) 
1 



9i±92 2 



Ti- 



ll! m — 2 In m In 2 + In p — 



+ -L In m g? ± g?g? - —In mg? ± g R g? + ' 



2tt 2 
1 



2tt 2 



2tt 2 

In 2 m — 2 lnm( In 2 + lnp ) 



2vr 2 
9?J R2 



In m — 2 In m In 2 + In p — 



4 



2% 2 



lnm(g^ + gf 



2g R2 



9?f 



(20) 



From Eqs. ( p0|) and fll8"|) we derive (3f]_(g R ). The final result is 

Px±(9 R ) = l ~9 R \9^ + ^(9 R \9?± + 9?f) + 0(g% (21) 

It can be noted that p does not appear in Eq. (|21|), as expected. A crucial use in deriving 
Eq. (j2lD is made of Eq. fll9|) , which is responsible for the cancellations expected from the 
exact solutions of the Luttinger and Mattis modelsi'SS and for the anomalous behavior 
of the theory. Of course the anomalous exponent rj = ji(g*) derived from Eq. ([19]) when 
< g\ <C 1 is in agreement with the exact solution of the Luttinger model, where g\ = 0, 

34 = 0. 

Equation (plD is the same as Eq. (P) for the g±± coupling. For g% and gn the same 
conclusion holds: the beta function of the GML method is recovered. The present CS 
approach, admittedly too involved, has perhaps the only value in that no use is made of 
approximate multiplicative relations. 



IV. THE WILSON APPROACH 



The multiscale formulation of the Wilson RGlil is particularly well suited to study the 
running of the coupling constants by discrete steps. The application of this method to 
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interacting one-dimensional fermionic systems started with Refs. |^,|ll| and was thoroughly 
developed and applied to various problems0HO. Here we give a short and simplified account 
of the method and refer to the cited papers for the details. 

In the coordinate space the free propagator G u (x) for u> particles (again uj = ±1 and 
Vf = 1) is: 

]_ r e -i{koxo+kixi) 

G u {x) = ~tt \o / dkodki — - — ■ — — . 

Actually it is not necessary to start with a kinetic term linearized around the Fermi surface: 
the RG can deal with realistic quadratic dispersion relations!. This simplification is however 
inessential for our purposes. Let p be an arbitrary momentum scale which for instance may 
be chosen equal to the inverse of the range of the potential. The propagator is decomposed 
in the sum 



G u {x) = £ Gl h \x) 



with 



G«(x)=^/d*e 



-ikx x ° 



(2tt) 2 J —iko + tok\ 



Gl h \x)= j^-Jdk 



-ikx 



X 



(27r) 2 J —ik + ooki 

e -p- 2 l- 2h {kl+kl) _ e -p- a 7- ak+2 (*g+^) 



(22) 



where h < 0, 7 > 1 and kx = koXo + k\Xx. This decomposition divides the u.v. from the 
i.r. singularity of the propagator: G^ is singular in the u.v. while Y4i=-oo Gffl = G l J' 
is singular in the i.r. It is important to note that Gffl(x), the propagator on scale h, for 
h < has an u.v. and an i.r. cutoff: GW(i) is essentially different from only for x ~ r )~ h 
{k ~ 7^ in momentum space). 

One imagines that this decomposition stems from a similar decomposition of the fields: 

tx,u),u / j Tx,ui,a 

h=—oo 

such that the pairings in the Grassmannian Wick rule are 

= Su>,u>>$<T,<T>8h,h'Gt( x - y)- 

We are interested to study the i.r. effective potential arising from the integration of the 
u.v. component ipffl from the effective potential V e s(ip) defined by: 

e~ v ^ =jfj P(#) e~ v ^\ 

where M is a normalization constant and V is the interaction potential. The ultraviolet 
integration was actually performed for the spinless mo del0. In the following we suppose to 
start directly with V^ ' . 
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The core of the method consists of a procedure that, integrating out the fields from 
the higher to the lower scales h (h — > — oo), constructs a well defined dynamical system of 
running coupling constants gh, whose iteration map is the beta functional. 

The operators C and 1Z = 1 — £ are introduced. 1Z is the usual renormalization operator 
of the BPHZ scheme: its action on a given vertex T, in momentum space for instance, is 
given by 1Z(F) = F — t r (T), where t r denotes the Taylor series with respect to the external 
momenta of F up to order D(F), if D(F) is the F superficial degree of divergence. Of course 

An /'('')• 

The couplings gh on a given scale h are defined by an inductive scheme. Let us assume we 
have constructed the effective potential V^(tp^ h \ g h+1 , . . . ,g ) on scale h, where ipj^^ = 
J2 n <h an d 9h+i, ■ ■ ■ , go are the previously defined couplings on scales h + 1, . . . , 0, . We 
define 

V ih \^ h \g h )=CV^{^ h \g h+l ,...,g»). 

The previous relation introduces the gh and relates them to the gh+i, ■ ■ - ,9o through the 
beta functional B h : g h = g h+1 + B h (g h+1 , . . . ,g ). The effective potential V^ -1 ) on scale 
/i — 1 is defined by 

e _y(h-i)(^,(<h-i)) _ 

J_ j p( d ^(h)^ e -CVW(^ h) )-KVW(^ h) )_ (23) 

Of course V^ 1 ^ = V^ 1 ' (ip^" 1 ' , gh, ■ ■ ■ , go)- The procedure is then iterated. The starting 
point is given by the couplings go of CV^ . The final goal is to find a region in the space of 
parameters g where each initial value generates a trajectory g h = g h+ i + B h (g h+1 , . . . ,g ) 
such that the Schwinger functions are analytic in the g%. 

Unfortunately this scheme in our problem requires emendation. From the second order 
result it becomes clear that ah and (h grow too fast independently on the initial conditions. 
The point is that we know that the interacting propagator has an anomalous behavior: 
asymptotically for large distances it decays faster than the free propagator. The wave- 
function renormalization necessary to cure this problem is accomplished by an inductive 
procedure that redefines step by step the free measure of the functional integral and the 
couplings by the means of a sequence of parameters Zh with h = 0, — 1, . . .. Let us assume 
we have introduced Zh, Zh+i, ■ ■ ■ , Z and applied our procedure integrating out the scales 
from to h + 1 (h < 0). We get an effective potential V^- h ^ (different from V^- h \ defined 
by Eq. (f2|)). We denote with P Zh (ipi h) ), ^ h (^- h_1) ) and P Zh (^ h) ) the free measures 
with propagators, respectively, Gffl jZ h , G^ h ~^ /Z h and G^' jZ h , where the last one is the 
modified propagator on scale h and G^ h ~^ = J2i<h-i 

Gf«. V(h-l) jg defined by 

J PzM^-^PzM^) e-^^ h) \ (24) 

y( h -l)^ Z l/2^ { <h-l^ hag ^ form . 

v-(<fc-i)( Z i/^(<A~U) = (c + n^^iz'fy^) 
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d 2 k (</,_!) + (</!_!)_ 

2 ^fc,w,cr rfc,u;,cr 



(2vr 



/" ^ , (<h-l)+/ •; \ i(<h-l)- 

<h-i 2^ J (-^o) 



(2 



Now we add and subtract from V^ h 1 - ) the term oc ZhCh-ii't^a { u ^i)'4 , hua ^ an d insert 
the term oc ZhCh-i^t^a + w&i) ipku^a i n the free measure. Let PzA^uf 1 ^ 1 ^) be 

the measure changed this way. We define Zh-i = Zh(l + Ch-i) and write: 

Pz h {d^ h ~ l) ) e-^ h - 1] ^^ h ~^ = J P' Zk (d^ h -V) e- y ' {h - l) ^£ h ~^ (25) 

^.(^^-^^-.(^r 1 ^-^" 15 ^^^ (26) 



Pz^Wg^Pz^Wr*) e-^V^^- - \ (27) 



In Eq. (P^D yC 1-1 ) is obtained from V"^ -1 ) dropping the Ch-i term and substituting a^-i — 
(h-i for ah-i- In Eq. (|26 ) Pz^Adip^ 1 ') is the free measure with propagator G^ 1 ' / Zh-i 
defined such that the remaining part of the free measure is exactly P Zh _ 1 (#^ h " 2 )). Finally 
Eq. ( ^7|) defines V^ -1 ) and the r.h.s has the same structure of the r.h.s. of Eq. ( p4|) so 
the procedure may be reiterated. The starting point are Z = 1 and Gffl = G^\ The 
relations between the couplings of CV^- h ~ lS) (ip^ h ~^) , g[ h _i (i = 1,2,4), cth-i and (h-i, and 
the couplings of g^h-i and S^-i, are easily found: 

h-i = ^—(cLh-i - Zh-i) 

9i > h - x =\~^~) 9 ' l ' h - 1 z = 1 '2,4. 

Of course (h-i is no more present in V r (^ ft —i) (^p 1-1 )) , but Z^-x is introduced. The replace- 
ment — > 5h drastically improves the convergence properties in the limit h — > — oo. 

Now we have the full recipe to proceed. Needless to say a crucial use of the linked cluster 
theorem will be made. For brevity only the calculation for the g\ coupling is sketched (we 
take g x ± = g x \\ = g x ). 

Let Cij denote the loop of figure § 

Cy = = / (02 ^)G^(fc), (28) 

where i,j < and is the propagator in momentum space on scale i. It is immediate to 
verify that the r.h.s. of ( ^8|) does not depend on uj nor on the scale p (see (0)) and that it 
is a function only of the difference i — j. In particular = C , z < 0. The second order 
calculation gives for g^: 
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9i,h-i = 9i,h + -K{c Kh g\ h + 2 E Chj g\ 

V h<j<0 



(29) 



K is a combinatorial factor: = 4. The previous relations give the second order approx- 
imation for the beta functional. Since we are interested in the beta function, in Eq. fl2"9"|) 
we write the g±j for j > h as functions of g^. This inversion generates a correction to the 
higher orders, in particular to the third order. We have (up to O(gf) terms): 



gi,h-i = gi,h + 7; K ( c h,h + 2 E c h,j)gl,h 

h<j<0 ' 

-k 2 e cj E + 2 E X o,A'/ : ,v 



(30) 



h<i<0 v h<j<i 

Using Eq. @ for C?(J0 we find: 



h<k<i k<j<0 



C h)h + 2 E C fc)i = ^-[+ln7-ln(l + 
+ ln(l+ 7 - 2 ^ 1 1 



,2-2h 



7 



27T 



ln 7 + 0( 7 2ft ). 



(31) 



From Eqs. ( pCf ) and ( |3~1~D we get the second order discrete beta function, up to r y 2h terms 
(let's remember that 7 > 1 and that we are interested in the h — > —00 limit): 



9l,h-l — 9l,h 



l,11 -9l h + 0(9 3 ). 



71 



The second line of Eq. (|30|) gives the corrections to the third order result. We find: 



E C h,i ( E C 3,j + 2 E X ^Vi 

h<i<0 ^ h<j<i h<k<i k<j<0 

1 



2tt 



/i<i<0 

where the r.h.s. is easily calculated: 



E In (l+7 24 ) + (*-/>) In 
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E ^ln(l + 7 2i ) < V E 1 



h<i<0 

/i<i<0 ^" 



h<i<0 



(32) 
(33) 



A is a constant. We neglect the r.h.s of Eq. (|32|). Equation ( p3|) too is derived neglecting 
terms exponentially vanishing with h. Putting all together we find a correction c to the 
third order given by: 

2hi2hi7 , 
r = 2 9i, h - 



7T* 



The computation of the third order is tedious. We will limit ourselves to note that exists 
a contribution oc In 2 In 7 of two-loop diagrams that cancels exactly c. This is important 
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because there is no such term in Eqs. (|6|), (0) or fl2ip . It comes from the diagrams D\ and 
D 2 of figure |3], which are related by D\ = —2D 2 . We will consider the simpler D\. The 
contributions to D\ given by 

y Ch,h^Ch^ + y y 20^20^ 

h<i<0 h<i<0 h<j<i 

+ /"I X/ Ch,i2Cj,k + /J y 2Chj2Ck,i 

h<i<0h<j,k<l h<i<0 h<j<i h<k<j 

amount to 

_i_ ln21n7 __L^ 2Cjikln(l + f .). (34) 

The second term of Eq. ( p4| ) can be neglected as for Eq. (p2|) and the first one gives the 
desired cancellation (taking into account the combinatorial factors). The final result is the 
same as Eq. (0) or (|2l| ) with g 1{ = g x ± = g 1 . Again we find the fixed point g* of Eq. (||) 
and we recover Eqs. (^]) in the limit 7 — > 1. 



V. CONCLUSIONS 

We have computed the third order (two-loop) approximation of the beta function for 
a one dimensional model of interacting fermions, aiming in particular to study the case of 
attractive interaction. An existing result!!, derived making use of tacitly assumed approx- 
imations, pointed out a 0(1) fixed point. We tried to support this conclusion setting a 
Callan-Symanzik approach and using the Wilson RG formulated as in Refs. [29|,[^,|ll]. In 
each case we recovered the aforementioned result. 

An attempt to pursue further the examination of the problem was made in Refs. |20|-f22 
where the fourth-order approximation, which is by no means universal, was computed. We 
propose a different approach focused on the study of the dependence on 7, the rescaling 
factor of the RG group. A similar idea was discussed in Ref. [32], where the dependence of 
the fixed points on the parameters of the RG was analyzed. Our simple idea is that the 
stability of the result with respect to 7 should indicate how reliable one should consider 
the perturbative result. It was expected a third order result dependent on 7 but we found 
a too strong dependence: the fixed point happens to change sign if 7 > yfe, which is still 
~ 1 (of course taking 7 > 1 and, at the same time, truncating at the third order would be 
questionable). Nothing similar happens to the nontrivial fixed points for repulsive interaction 
{g\ = 0, g 2 }, which are in some sense insensitive to the value of 7. 

Which conclusions can be drawn? It is useful to compare the one dimensional interacting 
Fermi Gas with the well known Kondo problem. It can be noted that the scaling equation 
for g\ (Fermi Gas) and the one for the impurity coupling A (Kondo model) have the same 



structure (see e.g. Ref. [33] for a review in the modern language of Conformal Field Theory). 
The Kondo effect was thoroughly investigated. The ferromagnetic case corresponds to the 
Fermi gas with repulsive interaction (gi > 0): the RG flow is such that A — > 0. If the 
coupling is antiferromagnetic the system flows toward a strong coupling phase (A — > 00)0. 
Moreover the infrared divergences induce a scale, the Kondo temperature Tk, characterizing 
the low energy physics. 
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In our case the particular instability of the perturbative result with respect to 7, a 
dimensionless parameter without a physical meaning, should indicate that the RG flow does 
not actually stop at a finite value and suggests a conclusion similar to the previous one. In 
our case the characteristic scale should be a gap A for the spin degrees of freedom, whose 
expression, according to Ref. should have for small coupling an expression of the type 
A oc y/giexp(-l/gi). 
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FIGURES 

FIG. 1. One loop diagram contributing to the four-point vertex function. The full (dashed) 
lines represent right (left) moving fermions. The value of the graph is proportional to 

-2^(10 + ^(1 + !). 

FIG. 2. One-loop contribution defining Q j . k is the internal momentum. The two propagators 
are on scales % and j. 

FIG. 3. Two-loop graphs contributing to the beta function: D\ diagrams (a) and D2 (b) give 
simply related contributes: D\ = — 2D2. 
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